Fe clusters on Ni and Cu: size and shape dependence of the spin moment 
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We present ab-initio calculations of the electronic structure of small Fe clusters (1-9 atoms) on 
Ni(001), Ni(lll), Cu(001) and Cu(lll) surfaces. Emphasis is given on the spin moments, and 
their dependence on cluster size and shape. We derive a simple quantitative rule, which relates the 
moment of each Fe atom linearly to its coordination number. Thus, for an arbitrary Fe cluster the 
spin moment of the cluster and of the individual Fe atoms can be readily found if the positions of 
the atoms are known. 



Small atomic clusters on surfaces constitute an ex- 
tremely interesting subject, as their electronic structure 
lies in the transition from the behaviour in bulk to that 
of free molecules or atoms. In particular, the electronic 
structure is characteristic of the cluster in that it de- 
pends on the cluster atoms, shape, size, and orientation, 
as well as on the substrate material on which the cluster 
is deposited. Thus it is terra incognita for the electronic 
properties. Moreover, since they are not hidden in the 
bulk, such structures are directly accessible by a number 
of experimental techniques, such as scanning tunneling 
microscopy, x-ray magnetic circular dichroism, photoe- 
mission, etc. 

In the area of magnetism, small clusters of magnetic 
atoms on surfaces are expected to experience an enhance- 
ment of the magnetic properties. As the cluster size de- 
creases and the average coordination of the atoms be- 
comes smaller, the decreased hybridisation of the atomic 
wavefunctions should lead to more pronounced mag- 
netic effects. Such trends with cluster size have been 
monitored in recent experimental works; for example, 
using x-ray magnetic circular dichroism, Lau and co- 
workers |lj have focused on Fe/Ni(00T) systems, while 
Gambardella and collaborators Q have examined Co 
clusters on Pt(lll); and nuclear resonant scattering of 
synchrotron radiation has been employed for the study of 
Fe islands on W(110) ||. Also other techniques, such as 
perturbed angular correlation 4] or spin-polarized scan- 
ning tunneling spectroscopy |5|, allow the analysis of 
magnetic atoms on surfaces in the sub-monolayer regime, 
giving further impetus to the held. 

Similar to the ferromagnetic Fe, Co, and Ni clusters, 
V, Cr, and Mn clusters should also have a strong intra- 
atomic exchange field; but, in contrast to the ferromag- 
netic clusters, their inter-atomic exchange should be an- 
tiferromagnetic, and due to competing interactions the 
resulting magnetic order is in general non-collinear and 
depends on individual details. Thus it will differ among 
clusters, and to investigate the individual magnetic order 
of such clusters is a challenge for both experiment and 
theory. 

But even for homo-atomic and mass-selected clusters 
it is very difficult to address experimentally the magnetic 
properties of each one individually, let alone each atom in 



a cluster. In many cases just the mean moment of clusters 
of a particular size or the average mean magnetic moment 
per atom is found, averaged over an ensemble of clusters 
of the same size but different shapes. This motivated us 
to investigate the local and global spin moments of small 
Fe clusters (up to 9 atoms) of monoatomic height and 
varying shape and size with a material specific theory 
of predictive power, namely density functional theory. 
In this paper we summarise our results by establishing 
a simple quantitative rule which relates the local spin 
moment of an Fe atom in the cluster linearly to the co- 
ordination number of its nearest neighbour atoms, with 
little influence of the structure of the atoms in the rest 
of the cluster. In this way for any arbitrary Fe cluster 
the spin moment of the cluster and of the individual Fe 
atoms can be readily found if the positions of the atoms 
are known. Even if the cluster shapes and sizes are not 
exactly known experimentally, by making a reasonable 
assumption of the variation of sizes and shapes the pro- 
posed rule allows an estimation of the variation of the 
total moment across the deposited clusters. Motivated 
by the work of Lau et al. l| we concentrated our at- 
tention on Fe clusters on the (001) and (111) surfaces 
of the ferromagnetic substrate Ni and the non-magnetic 
substrate Cu. 

First-principles calculations of clusters on surfaces can 
be computationally quite demanding, considering that 
such systems break the translational symmetry in all 
three directions. Most ab-initio methods explicitly ex- 
ploit the translational symmetry and the clusters on sur- 
faces are approximated by huge, computationally expen- 
sive supercells which are periodically repeated. So far 
they have been applied mostly for magnetic adatoms and 
chains |(|. A breakthrough in the treatment of mag- 
netic clusters was the development of the Green function 
method of Korringa, Kohn and Rostoker 0, for the 
embedding of clusters at surfaces. Subsequent applica- 
tions 0, fiol ITU ] established this method as a powerful 
tool. 

Our calculations are based on density-functional the- 
ory in the local spin-density approximation 

(LSDA). 

The Green function method of Korringa, Kohn, and Ros- 
toker (KKR) is employed to determine the spin den- 
sity and effective potential. Within the method, the 
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Green function of the perturbed system, i.e., surface plus 
cluster, is related to the one of the "host" reference sys- 
tem (the clean substrate) by a Dyson equation. In this 
way the correct host boundary conditions are included 
automatically in the Green function, thus no supercell 
construction is needed. In order to allow for a screening 
of the perturbation induced by the cluster, at least the 
first neighbouring sites of all cluster atoms were consid- 
ered in the self-consistency cycle. The angular momen- 
tum expansion of the Green functions was truncated at 
( max = 3. Tests have shown that these parameters give 
reliable results for these systems. 

A full-potential approach with the correct description 
of the atomic cells |l4j was used. Since we are focus- 
ing here on trend calculations, the atoms in the clusters 
were situated in the unrelaxed lattice positions, using the 
LSDA equilibrium lattice parameter of 6.46 au (3.417 A) 
for Ni and 6.63 au (3.507 A) for Cu [H. 

The clusters considered on the Ni(001) surface are 
shown schematically in Figure^viewed from the top (all 
atoms lie on the surface) . The view is adapted to surface 
geometry, meaning that it is rotated by 45° with respect 
to the in-plane fee cubic axes. The smallest cluster is a 
single Fe adatom, while the largest consists of 9 Fe atoms. 
In each atom, the calculated spin moment is written, and 
the average (per atom) moment of each cluster is also 
given. The Fe moment is always ferromagnetically cou- 
pled to the Ni substrate moment and to the Fe moments 
within the clusters. Already at a first glance it is obvious 
that the average moment of the clusters depends on the 
cluster size. The single adatom has manifestly the high- 
est moment (3.24 /is), while the 9- atom cluster shows a 
lower average moment of 2.85 \xb- 

This behaviour is expected on the grounds of hybridisa- 
tion of the atomic d levels with the neighbours. Atoms in 
larger clusters have, on the average, higher coordination, 
thus their d wavefunctions are more hybridised; this leads 
to lesser localisation and lesser tendency to magnetism. 

To pursue this idea further, we tried to correlate the 
local atomic spin moment to the nearest neighbour co- 
ordination of each atom, irrespective of the form or size 
of the cluster. For instance, let us focus on all Fe atoms 
which have only one first Fe neighbour, i.e., N c = 1 (the 
coordination to the substrate is the same, N s = 4, for all 
Fe atoms). Such atoms appear in the clusters with size 
2, 3, 4, 5, and 7; there are, in total, 10 such examples 
(having excluded cases which are trivially equivalent by 
symmetry). All of them have spin moments ranging in 
the small interval between 3.10 and 3.13 hb- Similarly, 
for the Fe atoms with two Fe neighbours the spin mo- 
ment ranges from 2.91 to 2.97 [is- Collecting all possible 
cases, from N c = (single adatom) to N c = 4, we present 
the results in Figure^. 

A surprising feature is the almost linear dependence 
of the spin moment on the coordination number. In the 
general case, for arbitrary magnetic atoms, this is not 
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FIG. 1: A: Spin moment (in [ib) of atoms in 14 different Fe 
clusters ranging from 1 to 9 atoms on the Ni(001) surface, 
and average (per atom) moment of the clusters (the view is 
surface-adapted, i.e., rotated by 45° with respect to the in- 
plane fee cubic axes; the clusters are viewed from the top, i.e., 
all atoms lie on the surface). B: Linear trend for the atomic 
Fe spin moment as function of the coordination number of 
nearest Fe neighbours. 



expected. For example, Ref.0 shows a study of the mag- 
netism of small 4d atom clusters on Ag(001). There, 
it is reported that the susceptibility is highly non-local, 
resulting even in an increase of the spin moment with 
increasing coordination. This is related to the larger ex- 
tent of the 4d wavefunctions compared to the 3d ones of 
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Fe. In our case, Fe has a strong intra-atomic exchange 
field, arising from rather localised 3d wavefunctions and 
their positioning with respect to the Fermi level Ep. In 
a non-magnetic picture, Ep is well within the d virtual 
bound state, and the non-magnetic density of states at 
Ep, n(Ep), is always well above the transition point to 
the magnetic state: the Stoner criterion / • n(Ep) > 1 
(with / being the exchange integral) is safely fulfilled. 
Thus, the hybridisation of the 3d levels with the neigh- 
bours' wavefunctions does not affect the nonmagnetic 
n(Ep) critically. With increasing coordination, n(Ep) 
is gradually lowered, and so is the intra-atomic exchange 
field. The tiny variation of the local moments for Fe 
atoms with the same nearest neighbour coordination but 
being in different clusters may arise from intra-cluster 
interference effects or indirectly via the different polari- 
sation of the substrate. 

Preliminary results on Co clusters, for which the 3d 
virtual bound state is lower in energy, show that this 
linear behaviour holds to a lesser extent, and much less 
so for Ni, where the intra-atomic exchange field is even 
weaker [l6|. A similar difference among Fe, Co, and Ni 
clusters on Ag(OOl) has been shown in Ref. an d 0- 

On the other hand, Mn and Cr should have a strong 
intra-atomic exchange field, but they present tendency 
for antiferromagnetism or even non-collinear magnetic 
structures 16]. In this respect, Fe is expected to be per- 
haps unique among the transition elements in showing 
such a clear cut linear trend. 

Similar is the case of Fe clusters on the more compact 
Ni(lll) substrate. The considered clusters, sized up to 
8 atoms, are shown together with the spin moments in 
Figure|2K. Since the Fe adatom on the (111) surface has 
only 3 Ni neighbours instead of 4 on the (001) surface, 
its spin moment is higher than on the (001) surface. On 
the other hand, the number of Fe neighbours (in-plane) 
can increase up to N c — 6, which makes the Fe moment 
decrease up to 2.45 /is; the collected statistics (spin mo- 
ments vs. N c ) are presented in Figure 03. Again we see 
a linear behaviour of the moment vs. N c . 

Next we turn to Fe clusters on Cu(001) and Cu(lll) 
surfaces, for which we considered the same types of clus- 
ters as in the case of Ni. Since the Cu lattice parameter 
is slightly larger than the one of Ni, one might expect 
a lesser degree of hybridisation and therefore increased 
spin moments at the Fe atoms. Nevertheless, the mo- 
ments are slightly lower, presumably because in the case 
of Ni the magnetisation of the substrate assists the spin 
polarisation of the cluster. The results for clusters on 
Cu(001) and on Cu(lll) are shown in Figure|3J\ an( i 
respectively. The result for the single adatom, showing 
again the highest moment, should not be taken literally: 
it is well known that magnetic impurity atoms in the 
bulk and on the surface of nonmagnetic hosts enter the 
Kondo regime at low temperatures, whence the spin mo- 
ment shows quantum fluctuations and has a net average 



(A) Moments in Fe clusters at Ni(111) surface 
3A 3B 3C 

3.00 i 3.00 ) /-~V^/-^\ (2.95M.17S 

3.17) 

Average=3.16 Average=3.10 
4A ; 3 .i 9 ) " B SA i3 .04i 5B 

2.89 ! 3.04 ; : 2 .8^) ^2-74; 3.04; 
3-05! 2.89: 3.03 5 , ^ 

Average=2.97 Average=3.01 Average=2.92 Average=2.96 

(3.04) 6C 

I2.89A2.78 







3.03J 
Average=2.83 




1 2 3 4 5 
Number of Fe neighbours 

FIG. 2: A: Spin moments (in /is) of atoms in 15 different 
Fe clusters ranging from 1 to 8 atoms on the Ni(lll) surface, 
and average (per atom) moments of the clusters (the clusters 
are viewed from the top, i.e., all atoms lie on the surface). 
B: Linear trend for the atomic Fe spin moment as function of 
the coordination number of Fe neighbours. 



of zero. This regime is not accessible by a mean field 
theory as the LSDA. Thus the result shown here for the 
adatom should be considered either as the trend of the 
LSDA calculation, or what one expects above the Kondo 
temperature. For magnetic clusters (dimers and larger) 
the quantum fluctuations of the moment are strongly 
suppressed and the Kondo temperature decreases drasti- 
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FIG. 3: Linear trend for the atomic Fe spin moment as func- 
tion of the coordination in Fe clusters on Cu surfaces. A: on 
Cu(001), and B: on Cu(lll). 



cally with cluster size; thus the results given here for the 
clusters are realistic. 

The linear dependence of the spin moments is again 
evident, except perhaps the result for the single adatom 
(N c = 0). A fit of the calculated data, using a linear 
equation of the form 

M = -aN c + b, (1) 

where M is the local spin moment, results in values for 
a and b given in Tabled We note that calculations on 
Fe nanostructures on Ag(001), reported for instance in 
Refs. and [12 (using the KKR method), and [13 (us- 
ing the SIESTA code) , show a similar trend, although in 
these works fewer kinds of clusters are presented. 

In summary, we have presented first-principles calcu- 
lations of small Fe clusters on Ni and Cu (001) and (111) 
surfaces. The average (per atom) spin moments of the 
clusters are reduced with cluster size, mainly because the 
higher average coordination of the Fe atoms increases the 
hybridisation of the d states. The important parameter 
turns out to be the local coordination of each distinct Fe 
atom, irrespectively of the exact size or shape of the clus- 
ter. A linear relation of the coordination to the atomic 
spin moment has been found in a good approximation. 



Surf, ci(hb) Mjib) 
Ni(001) 0.17 3.29 
Ni(lll) 0.15 3.33 
Cu(001) 0.15 3.21 
Cu(lll) 0.12 3.23 

TABLE I: Parameters a and b from the linear fit of the spin 
moment according to Eq. (1). 
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